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INDIVIDUAL VERSUS CLUSTER RECOVERIES WITHIN A 
SPATIALLY STRUCTURED POPULATION 

By L. Belhadji and N. Lanchier 

CNRS — Universite de Rouen 

Stochastic modeling of disease dynamics has had a long tradition. 
Among the first epidemic models including a spatial structure in the 
form of local interactions is the contact process. In this article we 
investigate two extensions of the contact process describing the course 
of a single disease within a spatially structured human population 
distributed in social clusters. That is, each site of the c?-dimensional 
integer lattice is occupied by a cluster of individuals; each individual 
can be healthy or infected. The evolution of the disease depends on 
three parameters, namely the outside infection rate which models 
the interactions between the clusters, the within infection rate which 
takes into account the repeated contacts between individuals in the 
same cluster, and the size of each social cluster. For the first model, we 
assume cluster recoveries, while individual recoveries are assumed for 
the second one. The aim is to investigate the existence of nontrivial 
stationary distributions for both processes depending on the value of 
each of the three parameters. Our results show that the probability 
of an epidemic strongly depends on the recovery mechanism. 

1. Introduction. To study the course of a disease within a spatially struc- 
tured population, Harris [8] introduced a model known as the basic contact 
process. Each site of the d-dimensional integer lattice is occupied by an indi- 
vidual either healthy or infected; each individual gets infected at a rate that 
depends on the number of infected individuals in some interaction neighbor- 
hood. Including spatial structure in the form of local interactions has shown 
that, for a disease to spread, the infection rate needs to exceed a threshold 
that is greater than the one for a nonspatial population. The smaller the size 
of the interaction neighborhood, the greater the threshold for the disease to 
spread. The reason for this is the lack of sufficient numbers of susceptible 



Received June 2005; revised July 2005. 

AMS 2000 subject classifications. Primary 60K35, 82C22; secondary 92D30. 
Key words and phrases. Contact process, epidemic models, Richardson's model, ran- 
dom walk, rescaling argument, individual recovery, cluster recovery. 

This is an electronic reprint of the original article published by the 
Institute of Mathematical Statistics in The Annals of Applied Probability^ 
2006, Vol. 16, No. 1, 403-422. This reprint differs from the original in pagination 
and typographic detail. 



1 



2 



L. BELHADJI AND N. LANCHIER 



individuals near the location of a disease outbreak once the disease starts 
spreading. 

The study of the evolution of diseases within spatially structured pop- 
ulations, including extensions of the contact process, is widespread in the 
particle system literature. The first process we investigate in this article has 
been introduced by Schinazi [12], and will be referred to as the cluster re- 
covery process (CRP). His model studies the spread of an infectious disease 
such as tuberculosis within a population grouped in social clusters, each 
of the clusters having the same size. The second model we investigate as- 
sumes another recovery mechanism and will be referred to as the individual 
recovery process (IRP). 

For both CRP and IRP, the dynamics depends on three parameters, 
namely the outside infection rate A (the rate at which an individual in- 
fects healthy individuals of other clusters), the within infection rate (j) (the 
rate at which an individual infects healthy individuals present in the same 
cluster) and the cluster size k (which is constant regardless of the location of 
the cluster on the lattice). The population is divided into social clusters, the 
individuals belonging to the same cluster having repeated contacts and the 
individuals belonging to neighboring clusters having casual contacts only, 
which suggests that the infection spreads out faster vertically within the 
clusters than horizontally between the clusters. In particular, it is assumed 
for both models that, once a cluster has at least one infected individual, 
infections within the cluster are much more likely than additional infections 
from the outside so we neglect the latter. The cluster size can be seen as 
the mean number of individuals having sustained contacts with a given in- 
dividual. Even if, in a more realistic setting, this parameter should fluctuate 
significantly depending on the social customs of each of the individuals, it 
is assumed for technical reasons that all the clusters have the same size. In 
particular, whereas the infection rates A and (p are parameters linked to the 
nature of the disease, the cluster size k depends on the social customs of 
the individuals. The only difference between the CRP and the IRP is the 
recovery mechanism. For the CRP, all the infected individuals in a given 
cluster simultaneously recover at rate 1 thanks to an antidote. This applies 
to places where there is a good tracking system of infectious diseases so that, 
once an infected individual is discovered, its social cluster is rapidly tracked 
down [12]. For the IRP, we deal with the other extreme case when at most 
one infected individual recovers at once; that is, the tracking system is not 
effective enough and the infection can spread within a given cluster before 
it is detected. In particular, the CRP and the IRP can be considered as 
spatial stochastic models for the transmission of infectious diseases in de- 
veloped and developing countries, respectively. As we will see further, this 
difference implies that, the cluster size being fixed, an epidemic may occur 
for the IRP provided the outside infection rate is strictly positive. In the 
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CRP, the whole population always recovers if the outside infection rate is 
smaller than some critical value depending on the cluster size. Moreover, we 
obtain for the IRP that, when the within infection rate is greater than 1, 
even if the outside infection rate is low, an epidemic is possible provided the 
cluster size is large enough. The CRP exhibits the opposite behavior in the 
sense that, the within infection rate being fixed, if the outside infection rate 
is low, there is no epidemic whatever the cluster size. 

To figure out the differences between individual and cluster recoveries, 
we now introduce the explicit dynamics of both processes. The IRP is a 
continuous-time Markov process in which the state at time t is a function 

'■ — > {0, 1, . . . , k}, with K denoting the common size of the clusters. The 
cluster at site x € is said to be healthy at time t > if Ct{x) = 0, and 
infected otherwise. More precisely, (,t{x) indicates the number of infected 
individuals present in the cluster at time t >0. To take into account the 
outside infections, we also introduce an interaction neighborhood. For any 
X, z € Z'^, X ~ z indicates that site z is one of the 2d nearest neighbors of 
site X. Then, the state of site x flips according to the following transition 
rates: 

— > 1 at rate A^^(z), 

1 — >■ i + 1 at rate icj), i = 1, 2, . . . , k — 1, 

— 1 at rate i, i = l,2,...,K. 

In other words, a healthy cluster at site x gets infected, that is, the state 
of X flips from to 1, at rate A times the number of infected individuals 
present in the neighboring clusters. In other respects, if there are i infected 
individuals in the cluster, i = 1, 2, . . . , At — 1, then each of these individuals 
infects healthy individuals in the cluster at rate (p. Finally, each infected 
individual recovers at rate 1 regardless of the number of infected individuals 
in its cluster. 

The CRP is a Markov process rjt : Z'^ — > {0, 1, . . . , k}, with r]t{x) denot- 
ing the number of infected individuals at site x at time t >0, and whose 
dynamics is obtained by replacing the transitions i— >i — 1, i = 1,2, . . . , k, 
above by the transitions 

i — > at rate 1, i = 1, 2, . . . , k. 

That is, all the infected individuals in a given cluster are now simultaneously 
replaced by healthy ones at rate 1, the infection mechanism being unchanged. 
For more details, see [12]. 

The graphical representation. An argument of Harris [7] assures us of the 
existence and uniqueness of our spatially explicit, stochastic models. The 
idea is to construct the processes from collections of independent Poisson 
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processes, which is referred to as Harris' graphical representation. For each x, 
z G Z*^ with X ~ z and z = 1, 2, . . . , k, we let {T^'^'* : n > 1} denote the arrival 
times of independent Poisson processes with rate A, and draw an arrow 
labeled with an i from site x to site z at time T^''^'* to indicate that an outside 
infection may occur. To take into account the within infections, we introduce, 
for X S Z"^ and i = l,2,...,«; — 1, a further collection of independent Poisson 
processes, denoted by {C/^'*:n > 1}; each of them has rate (j). We put the 
symbol at (x, C/^'*) to indicate that an infection from the inside may occur. 
Finally, for each x € Z'^ and i = 1, 2, . . . , k, we let {1^'* : n > 1} be the arrival 
times of independent rate-1 Poisson processes, and put a Xj at site x to 
indicate that a recovery may occur. 

Given initial configurations and r/Oi and the graphical representation 
introduced above, the processes can be constructed as follows. If there are 
at least i infected individuals at site x at time T^'^'* and the cluster at z 
is healthy, then the state of z flips from to 1 for both processes. In other 
respects, if there are j infected individuals, i < j < k — 1, at site x at time 
[/^'\ then one more individual gets infected in the cluster, that is, the state 
of X flips from j to j + 1, for both processes. Finally, if there are j infected 
individuals, 1 < j < k, at site x at time V^'*, then the state of x flips from j 
to J — 1 if and only if i < j for the IRP, while it flips from j to if and only 
if i = 1 for the GRP. In particular, have no effect on the 

GRP. 

A nice feature of the graphical representation is that it allows us to couple 
several processes starting from different initial configurations, which can be 
done by using the same collections of Poisson processes. See [5], page 119 
and [10], page 32. 

The mean-field model of the IRP. To figure out the properties of our 
spatial models, the first step is to investigate their deterministic nonspatial 
versions called the mean-field model; that is, we assume that all sites are 
independent and the system is spatially homogeneous. This then results in 
a system of ordinary differential equations for the densities of healthy and 
infected clusters. The reason for introducing this model is that the existence 
of locally stable fixed points for the mean-field model may be symptomatic 
of the existence of stationary measures for the original spatial model. See [6] 
for different possible studies of a model. Let Ui denote the density of clusters 
with i infected individuals, z = 0, 1, . . . , k. The mean-field model of the IRP 
is then described by the following coupled system of ordinary differential 
equations: 
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"dt 

dui 
~dt 

dt 



Xuo'^iui - (1 + <f))ui + 2u2, 



i=l 



{i - l)(/>Ui_i - i{l + 4>)ui + + l)iij+i, i = 2, 3, . . . , K - 1, 



To find the condition for the existence of a nontrivial equifibrium, we start 
by setting the right-hand side of each of the equations equal to 0. Prom the 
last equation, we obtain 



K-l 



By induction, we get 



i-l 



Ui 



i'Ui-i, 



i — 2, 3, . . . , K, 



from which it follows that 



Ui 



-ui, z = 2, 3, . . . , K. 



Reporting in the first equation then leads to 

J re— 1 

-— = ui - A 2^ liQUi = 0. 
1=0 





Fig. 1. Solution curves of the mean-field model of the IRP in the case n = 2. Left: A = 1 
and <j> = 0.5. Right: A = 0.5 and (j) = 2. 
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The previous equation tells us that the condition for an epidemic to occur 
is given by 

(1) X{l + c^ + --- + <P^-^)>l. 

When (1) does not hold, the boundary equilibrium uq = 1 and Uj = for 
i = 1, 2, . . . , K is the unique fixed point of the system of ordinary differential 
equations. When (1) is satisfied, there is a nontrivial equilibrium character- 
ized by 

Uq = — X < } > and Ui = ui , i = 2, 3, . . . , k. 

^ l^^ J 

First of all, note that, when A > 1, condition (1) holds for all the values of 
> and K > 1. In other respects, when A < 1 and At > 2, there is a nontrivial 
equilibrium provided (j) is sufficiently large. In the same way, when A < 1 and 
(f)>l, there is a nontrivial equilibrium if the cluster size k is large enough. 
Moreover, numerical simulations indicate that, when condition (1) holds, 
the nontrivial equilibrium is locally stable, which suggests the existence of 
a nontrivial stationary distribution for the corresponding spatial model. See 
Figure 1 for a picture. For a description of the mean-field model of the CRP, 
see [12], Section 3. 

The stochastic models. We now discuss the effects of each of the three 
parameters, namely the outside infection rate A, the cluster size k, and the 
within infection rate (p, on the probability of an epidemic for both models. 
We will prove that the behavior predicted by the mean-field model above 
holds as well for the IRP, and provide comparisons between both processes. 
Except for Theorem 6 whose proof is carried out in this paper, all the results 
regarding the CRP can be found in [12], so we now focus especially on the 
IRP, but the following definitions and basic properties hold for the CRP as 
well. 

We say that an epidemic may occur for the IRP with parameters (A, k, (p) 

if 

-P°(|6I > 1 for alH>0) >0, 

where P° denotes the law of the process starting from one infected individual 
at the origin, and where is the total number of infected individuals 
at time t > 0. The previous property is equivalent to the existence of a 
nontrivial stationary distribution, where nontrivial means ^ 6o, the "all 0" 
configuration. We say, on the contrary, that there is no epidemic if 6q is the 
unique stationary distribution. 

As in [12], one of the keys of our results is monotonicity. A crucial feature 
of the CRP that ensures monotonicity is the fact that the transition i ^ 
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occurs at a constant rate (actually, any transition rate which is decreasing 
in i would work). For the IRP, this follows from the fact that we have at 
most one recovery at any fixed time. The following monotonicity result can 
be checked by using a coupling argument if we think of the process as being 
generated by Harris' graphical representation introduced above. 

Lemma 1.1. The IRP is attractive and monotone with respect to A, k 
and (p. 

Here attractivity means that if ^o(^) — ^o(-^) rr E Z*^ at time 0, 

then ^1 and can be constructed in the same probability space in such a 
way that 

P^'2(Vx € (x) < Ct{x)) = 1 for any t > 0, 

with P^'^ denoting the law of the coupled process starting from (Co^Co)- 

In order to introduce our results, we first observe that, when k= 1, each 
of the clusters is only in one of the two states = healthy or 1 = infected, 
and the process reduces to a basic contact process with parameter A. In 
this case, there exists a critical value Ac € (0, oo) such that if A < Ac, then 
the process converges in distribution to the "all 0" configuration; otherwise, 
an epidemic may occur. See [2] or [10], Theorem 2.25. This, together with 
Lemma 1.1, implies that, when A > Ac, an epidemic may occur for any k>1 
and > 0. This case corresponds to the case A > 1 of the mean-field model. 

To figure out the intuition behind our first result, we start by removing the 
interactions between clusters by setting A = 0. This makes the IRP a system 
of independent random walks with absorbing state 0; each of them repre- 
sents the number of infected individuals in the associated cluster. Since each 
of these random walks returns to with probability 1, the process converges 
in distribution to the "all 0" configuration. By relying on a perturbation ar- 
gument, the result can be extended to the region A > small. More precisely, 
we have the following 

Theorem 1. For all k>1 and <f>>0, there is \c{K,(j)) G (0,oo) such 
that if X < Ac(k, (/>) there is no epidemic for the IRP, while if X> \c{k,4>) 
an epidemic may occur. 

Note that, due to the monotonicity with respect to k, we get Xc{K,(j)) < 
\c{l,4>) = Xc- The analogue of Theorem 1 for the CRP is given by 

Theorem 2 ([12]). For all k> I and (j)>0, if X < Ac/k there is no 
epidemic for the CRP, while if X> Xc an epidemic may occur. 
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The next step is to set A < Ac and k>2, and to discuss the probabihty of 
an epidemic depending on the value of the within infection rate (p. To begin 
with, when (j) = the IRP reduces to a basic contact process with parameter 
A. In particular, since A < Ac, there is no epidemic. The other extreme formal 
case (j) = oo corresponds to a Richardson model with parameter Ak so that an 
epidemic may occur provided A > 0. To see that the previous two conclusions 
still hold when (j) G (0,00), we will first rely on a rescaling argument to 
estimate the rate of convergence of = 0) in the two limiting cases 

<j) = and (j) = 00, respectively, where denotes the law of the IRP. These 
estimates will have to be good enough so that a perturbation argument can 
be applied. This, together with Lemma 1.1, implies that 

Theorem 3. For all k>2 and A G (0, Ac), there is 0c(A, k) € (0, 00) such 
that if (j)< 4>c{^-, k) there is no epidemic for the IRP, while if 4> > 4>c{^, n) 
an epidemic may occur. 

The CRP exhibits a quite different behavior since such a critical value 
exists if and only if the outside infection rate belongs to {Xc/ k, Xc), that is: 

Theorem 4 ([12]). For all k>2 and AG {Xc/k,Xc), there is (j)c{X,K) G 
(0,00) such that if 4> < (pciX,^) there is no epidemic for the CRP, while if 
(j)> (pdXjK) an epidemic may occur. 

The last step is to investigate the effects of the cluster size k on the 
probability of an epidemic for both models. For the mean-field model of the 
IRP, we have seen that, in the case A > and > 1, condition (1) holds 
provided k is sufficiently large. The assumption (j)>l is to make sure that 

the sequence 1 + (p-l h cl)'^~^ diverges as k — > 00. The analogous result for 

the spatial version is given by the following 

Theorem 5. For all (j) > I and X G (0,Ac), there is Kc{X,(j)) > 2 such 
that if K < Kc{X,(j)) there is no epidemic for the IRP, while if k > Kc{X,(j)) 
an epidemic may occur. 

To understand the assumption (p > 1, we observe that when the inter- 
actions between the clusters are removed, that is, when A = 0, the process 
becomes a system of independent random walks with absorbing state 0. 
When (j) > 1 and only in this case, these random walks have a drift to the 
right so that the infection in a given cluster can persist a very long time. 
As we will see further, the previous observation is the key for proving Theo- 
rem 5. Our last result tells us that the CRP exhibits the opposite behavior 
in the sense that, (j) being fixed, when A is too small, there is no epidemic 
whatever the cluster size. More precisely, we have the following 
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Fig. 2. Phase diagram of the stochastic models. The continuous curve refers to the IRP 
and the dashed one to the CRP. For the IRP, Theorems 1 and 3 imply that the straight 
line \ — IS an asymptote of the phase transition curve, and Theorem 5 that, when > 1, 
Xc{K,(j)) — > as K —> oo. For the CRP, Theorem 6 implies that \c{i^,(f>) — > Xc{<i)) > Q as 
K — > 00. The monotonicity follows from Lemma 1.1 and Proposition 1 of [12] while we 
conjecture the continuity of the two curves. 

Theorem 6. For all (j)>0 there exists Xd^P) ^ (0, cjo) such that for any 
K > 1 there is no epidemic for the CRP provided A < Ac(</>). 

The key idea of the proof is that, contrary to the IRP in which the number 
of infected individuals in a given cluster performs a random walk with a 
drift to the right when <j) > 1, for the CRP, the clusters globally recover 
sufficiently often so that the number of infected individuals between two 
recoveries cannot exceed some threshold Ko{(j)) with probability close to 1. 
In particular, by taking A > smaller than some critical value depending 
on KQ{(j)), but independent of k, the disease dies out. All our results are 
summarized in Figure 2 above. 

The rest of the article is devoted to proofs. In Section 2 we rely on a 
rescaling argument to prove Theorem 1. In Section 3 we use the fact that, 
when A € (0, Ac), the contact process with parameter A converges to the 
"all 0" configuration while the Richardson model with growth rate Xk has 
a nontrivial stationary distribution, to deduce Theorem 3. The proof of 
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Theorem 5, partially based on random walks estimates, is carried out in 
Section 4. Finally, Section 5 is devoted to the proof of Theorem 6. 

2. Proof of Theorem 1. In Section 1 we have seen that when A > Ac an 
epidemic may occur whatever the size of the clusters and the value of the 
within infection rate. In view of the monotonicity of the IRP with respect 
to the outside infection rate A, if in addition we prove that, k > 1 and 
(/) > being fixed, there exists Aq > such that, when A < Aq, the "all 0" 
configuration is the unique stationary distribution, Theorem 1 will follow. 

In order to prove the existence of such a Aq > 0, we use a rescaling argu- 
ment to compare our stochastic process with an oriented percolation process 
on Q = Z,'^ X Z+ . The rescaling argument has been invented by Bramson and 
Durrett [3] and is reviewed in [5]. To make the connection between the par- 
ticle system and an oriented percolation process, the basic idea is to turn 
the space-time of the process into a brick wall, and to associate each brick 
with a certain good event Ez^n, {z,n) £ Q. For any {z,n) G we consider 
the vertical segment 

(2) S,^n = {{z,t):nT<t<{n + l)T} 
as well as the space-time cylinder 

(3) C^^n = {{x, t) G Z'^ X M+ : X ~ z and nT < i < (n + 1)T} 

where T is a large integer to be fixed later. We say that a site {z^n) € ^ is 
good if the cluster at site z is healthy at time {n + \)T. We want the good 
event Ez^n to be measurable with respect to the graphical representation in 
Sz,n U Cz^n and such that on Ez^n site {z,n) is good regardless of the value 
of ^nT{z) or the configuration of the system in the cylinder Cz^n- We want 
this event to have a probability close to 1. 

To define Ez^n we start by setting A = to make our process a collection 
of independent random walks. Moreover, to be in the worst case, we assume 
that £,nT{z) = K. Let Tq be the first time S,t{z) hits after time riT. Since 
K < oo, the stopping time Tq is finite with probability 1. This, together with 
the Beppo-Levi theorem, implies that, for any e > 0, there is a large enough 
T, fixed from now on, such that 

(4) P^{T^ -nT<T)>l-e/2, 

where denotes the law of the particle system. To define the event Ez^n, we 
first require that Tg < (n + 1)T (which can be defined, in terms of Poisson 
processes, as a certain succession of x's and •'s at site z) so that on Ez^n 
site {z,n) is good provided A = 0. To take into account the local interactions 
between the clusters when A > 0, we complete the definition of Ez^n by 
requiring no arrow to come from the cylinder Cz,n to the segment Sz,n- 
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Finally, now that T is fixed so that (4) holds, we can find Aq > so small 
that 

is good) > Pi,{E,^n) >l-e if A < Aq. 

To complete the comparison, we now position oriented edges between sites 
in Q in order to obtain an oriented percolation model. For {x^n)^ (y,^) G 
we draw an edge from (x,n) to {y,m) if and only if a; ~ y and n = m, or 
X = y and m = n + 1. In that way, we may define a 1-dependent oriented 
percolation process on Q for which sites are open with probability 1 — e. By 
choosing e > sufficiently small, the probability of a path of length n of 
closed sites within this percolation process decreases exponentially fast with 
n (see (8.2) in [1]). Since the existence of an infected site at time (n + 1)T 
implies the existence of a path of length n of closed sites. Theorem 1 follows. 

3. Proof of Theorem 3. In this section we first prove that, given A € 
(0, Ac) and k>2, for the spatial model with parameters (A, k, (/>), there is 
no epidemic when (p = while an epidemic may occur when (j) = oo. Then, 
by using a perturbation argument, we will deduce from the behavior of the 
IRP in the two extreme cases that there exist (f>i> and (j)2 < oo such that 
if (p < (f)i there is no epidemic, while if (p > (f)2 an epidemic may occur. In 
view of the monotonicity of the model with respect to the parameter (p, this 
will imply the existence of a i;^c S [i;^i,i;^2] such that Theorem 3 holds. 

To prove the existence of (pi > small, we let L be a large integer to be 
fixed later and, for any site {z,n) G Q, introduce the space-time regions 

A, ^n = {Lz,Ln) + {[-2L,2Lf X [0,2L]} 

and 

B, ,n = {Lz,Ln) + {[-L,Lf x [L,2L]}. 

A site {z,n) € ^ is said to be good if (,t{x) = for any (x,t) € i?z,ni that 
is, there is no infected individual in Bz^n- When A < Ac and (p = 0, the IRP 
reduces to a subcritical contact process, which implies, by Section 5 of [12], 
the existence of an event Fz^n measurable with respect to the graphical 
representation of in Az^n such that on Fz^n site {z,n) is good. Moreover, 
for any e > there exists a large enough L, fixed from now on, such that 

Pg((z,n) isgood)>P5(F,,„)>l-e/2 if = 0. 

Now that L is fixed, we can find > small such that there is no occurrence 
of Poisson processes with rate (p inside Az^n with probability at least 1 — e/2 
provided (p <(pi- In conclusion, 

^^((z, n) is good) > 1 — e ii(p<(pi. 
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To complete the comparison with oriented percolation, we now position ori- 
ented edges between sites in Q as follows: For {x,n), {y,m) G ^ we draw an 
oriented edge from (x, n) to {y, m) if and only if n < m and A^^n H ^y,m 7^ 0- 
In that way, we define a 2-dependent oriented percolation process on Q for 
which sites are open with probability 1 — e. Since e > can be chosen ar- 
bitrarily small, we can conclude, as in Section 2, that there is no epidemic 
when (p <4>i. 

To get the existence of i?!)2 < oo, we now deal with the limiting case (j) = oo 
for which our process becomes a Richardson model '■ — > {0, k} with 
growth rate Xk (see [11]). If we think of the process as a set- valued process in 
which the state at time t is the set of infected sites, and start the evolution 
with one infected cluster at the origin, then the following shape theorem 
holds ([4], Theorem 6, Chapter 1): there exists a convex set A such that, for 
any 5 > 0, there are constants Ci < oo and 71 > such that 

(5) P^iil - 6)tA C 6 C (1 + 6)tA) >1-Ci exp(-7it). 
Let L and F denote two large integers to be fixed later, set 

H = {{z, n) ^1? : z + n IS even and n > 0}, 
and, for z G Z, introduce the spatial box 

(6) B, = Lzei + [-L,L]'', 

where ei denotes the first unit vector of the d-dimensional lattice. We say 
that a site (z, n) G W is occupied if each of the clusters in the spatial box Bz 
has K infected individuals at time nVL. Then, the shape theorem (5) implies 
that, for any e > 0, we can pick L and F sufficiently large so that the set of 
occupied sites for the Richardson model dominates the set of wet sites in 
a 1-dependent oriented percolation process with parameter 1 — e/2. 

To extend the result to i;^ < 00 sufficiently large, we require the following 
good event, denoted by Gz,n- For any x € B,-! U B^+i, we follow the line 
{x} X [riTL, (n -|- 1)FL] by going forward in time. Each time we encounter 
a X , we require the next • to appear before a new x or /t-arrow is encoun- 
tered. This good event assures us that the model with parameters (A, k, cp) 
exhibits the same behavior as the Richardson model introduced above in the 
spatial box B^-i U B^+i from time nTL to time (n -|- 1)FL. To estimate the 
probability of our good event, we first observe that, in view of the param- 
eters of the Poisson processes involved in our graphical representation, the 
probability of going through a • before any x 's or K-arrow is given by 

{k- l)4> + 2dX + K' 

Let M denote the number of x's contained in {B^-i U -B^+i} x [riTL, (n + 
1)FL]. In view of the properties of the exponential distribution, we get 

E^{M) = k{4L + 1)(2L + if-'^TL 
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since there are (4L + 1)(2L + l)'^"^ sites in i?z_i U B^+i- Finally, by decom- 
posing the event to be estimated according to whether M > 2E^{M) or not, 
large deviation estimates imply that 

^'5(G.,n)>l-C2exp(-72rL) 

-2K(4L + l)(2L + l)^"irLx - + 



l)(h + 2d\ + K 



for appropriate C2 < oo and 72 > 0. In particular, by setting cf) = exp(rL), 
and then choosing L sufficiently large, P^{Gz^n) can be made greater than 
1 — e/2. Putting things together, it follows that the set of occupied sites 
for the model with parameters (A, k, (f)) dominates the set of wet sites in a 
1-dependent oriented percolation process with parameter 1 — e. 

To construct our nontrivial stationary distribution, we start from the 
"all k" configuration, run the process to time 5, take the Cesaro average of 
the distribution at times < s < 5" and extract a convergent subsequence. 
By Liggett ([9], Proposition 1.8) the limit /x is a stationary measure. Since 
percolation occurs with positive probability when e > is small enough, 

liminf ^((z, n) is occupied) > limmlPMz^n) is wet) > 0, 

n^oo n— >oo 

where Pe denotes the law of the oriented percolation process with parameter 
1 — e, which implies that ^ concentrates on configurations with infinitely 
many infected clusters. 

4. Proof of Theorem 5. In this section we prove that if A G (0, Ac) and 

(/) > 1 there exists a critical value > 2 such that if k < Kc the IRP converges 
to the "all 0" configuration, while if k > Kc there is a nontrivial equilibrium. 
First of all, if k = 1 the value of (p is irrelevant and the process reduces to 
a basic contact process with parameter A < Ac so that the "all 0" config- 
uration is the unique stationary distribution. In particular, in view of the 
monotonicity of the IRP with respect to k, it suffices to prove the existence 
of a large kq such that the process with parameters (A, kq, (j)) has a nontrivial 
stationary distribution. This will imply Theorem 5. 

To prove the existence of such a hq , the strategy is to compare the particle 
system with a 1-dependent oriented percolation process on Ti. To rescale the 
particle system, we let = k'^, and say that a site {z,n) S is occupied if 
at any time t € [nT^, {n + l)?'^] there is at least one infected individual at 
site zei, where ei denotes the first unit vector of the d-dimensional lattice. 
Theorem 5 then follows from the following. 



Lemma 4.1. Let A > and (j) > 1. Then, for any e > 0, there exists a 
large enough k such that the set of occupied sites dominates the set of wet 
sites in a 1-dependent oriented percolation process with parameter 1 — e. 
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The first step is to investigate the behavior of an isoiated cluster when 
the outside infection rate is set to 0. More precisely, we need some estimates 
on the first recovery time of a large cluster, the main ingredient (given in 
Lemma 4.4) being that, with probability close to 1, the infection dies out 
quickly or persists at least 2Tk units of time. 

Random walk estimates. To figure out the behavior of an isolated cluster, 
we first consider the continuous-time random walk G {0,1,..., k} that 
makes transitions 

i + 1, at rate id>, , , ^ • ^ 

, . whenl<z<K — 1, 

I — I, at rate z, 

and K — 5- K — 1 at rate k. In particular, Xt is equal in distribution to ^t(x), 
the number of infected individuals at site x € Z'^, when A = 0. To estimate 
the recovery time of a given cluster, we also introduce the stopping time 

TO = inf {t >0:Xt = 0} 

and the good event 

= {Xt = K for some t < 2Tk}. 

The aim is to prove that, with probability close to 1 when k is large, tq > 2Ti^ 
on the event while the random walk returns quickly to on the event 
Q'f^. By convention, all through this section the subscripts and superscripts 
on the probabilities denote the process and its initial state respectively (e.g., 
for the law of the process Xt starting from Xo = k). 

Lemma 4.2. Let (/)> 1. Then there exist C3 < oo and 73 > such that 
P^{to < 2T^) < AK<l,T^r^ + C73exp(-73r«). 

Proof. The basic idea is to deduce the result from similar estimates 
for the asymmetric discrete-time random walk that makes transitions 

i + 1, with probability (/)/(l + (/>), < ^ _ 1 

1, with probability 1/(1 + (/>), wnen ^ ^ k i, 

and K — > K — 1 with probability 1. Let 1^, G Z be the backward random walk 
with transitions 

1 + 1, with probability + (p), 

2 — 1, with probability (/)/(l + 

and, for any i < 0, pi the probability that 1^ = for some n > when 
starting from Yq = i. Decomposing according to whether Yn first jumps to 
i — 1 or i + 1 leads to 

1+0 1+0 
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Since — > as i — > — oo, it follows that pi = (p'^, which implies that, for any 
n>0, 

Pxi^n < 0) < Py'^iYj = for some i > 0) = = 

Let Kt denote the number of times the process jumps by time t. Since 
makes transitions at rate at most K.(j), large deviation estimates imply that 

P^iKt > 2K(Pt) < C3exp(-73t/2) 

for appropriate C3 < oo and 73 > 0. Putting things together and decompos- 
ing according to whether K2Tk is smaller or greater than 4:K(pTi^, we obtain 

< PxiXn < for some n < Ak^T^) + Pj^(i^2T„ > 4k(/.T^) 

< ^ Pl(X„<0) + C3exp(-73r«) 

n=l 

< 4K0r«,^-'^ + C3 exp(-73T^). 
This completes the proof of the lemma. □ 

Lemma 4.3. Let (p > 1. Then there exist C4 < 00 and 74 > such that 
for any t < IT^ 

Pi.(ro>t;0^)<C4exp(-74t). 

Proof. To begin with, we prove the analogous result for the asymmetric 
random walk Zn (zT^ that makes transitions 

i + I, with probability 0/(1+0), 
i — 1, with probability 1/(1 + 0). 

Let o"o denote the first time Zn = and p = 0/(1 + 0). First of all, since 

E^z[Zn-Zo] = {2p-l)n, 

large deviation estimates on Zn imply that, for any 5 > 0, 

Pl{Zn <{2p-l- 5)n) < C5 exp(-75n) 

for appropriate C5 < 00 and 75 > 0. In other respects, since > 1, and so 
p > 0.5, we can find a sufficiently small 6 > such that 2p — 1 — 6 > 0, which 
implies that 

00 

Pl(n<(Jo<oo)= ^1(^0 = ^) 

k=n+l 
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oo 

< Yl PziZk<0) 

k=n+l 

oo 

<C'5 ^ exp(-75A:) 

k=n+l 

< C6exp(-75n). 

To deduce the result for Xg, we first observe that, on and by time 
min(To, 2rK), the continuous-time random walk is just a time change 
of Zn. Since jumps at rate at least 1 by time tq, the probability that 
Kf <t (where Kt denotes the number of jumps by time t) when tq > t is 
smaller than C7exp(— 77t). In conclusion, for any t < 2Tk, we get 

Pi(ro > t; ni) < Pi,(ro > t; ni;Kt >t)+ P),{ni; Kt < t) 

< < ao < oo) + Pkini; Kt < t) 

< C6exp(-75t) + C7exp(-77t) < C4exp(-74t) 
for appropriate C4 < cx3 and 74 > 0. This completes the proof. □ 

Lemma 4.4. Let (j)> 1 and = k^. Then there exist Cg < 00 and 78 > 
such that 

P^{t < To < 2T«) < C8exp(-78K) + C8exp(-780. 

Proof. To begin with, we observe that, on the event Or, the first time 
hits when ^0 = is bounded in distribution by the first time it hits 
when = 1 so that 

Pi (to < 2r,; < Pj^(To < 2T,). 

In particular, by decomposing the event to be estimated according to whether 
il^ occurs or not, it follows from Lemma 4.2 and Lemma 4.3 that 

Pi(t < TO < 2r,) < Pj^it < TO < 2T,; n^) + P^{t < TO < 2T^; n^) 

< PkiTo < n,) + Pi.(To > t; Ql) 

< 4k^(I)^~'^ + C3 exp(-73K2) + C4 exp(-74t) 

< Cs exp(-78K) + Cg exp(-78t) 

for suitable Cs < 00 and 73 > 0. This completes the proof of the lemma. □ 

Proof of Lemma 4.1. To deduce Lemma 4.1, we show that if the 
cluster at site has at least one infected individual from time to time T^, 
that is, site (0, 0) is occupied, then this individual gives birth at site 2; ~ 
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and by time to an infection that survives at least 2T^ units of time, so 
that (z, 1) is occupied. The probabiUty of this event has to be greater than 
1 — e for K sufficiently large. 

Each time the infected individual at site gives rise to a new infection at 
site z, we call this infection a strong infection if it lives at least 2T^ units 
of time (at site z). Let denote the first time a strong infection is born at 
site z, that is, 

fj^ = inf{t > : Cs(-z) / for any t < s < t + 2T^}. 

Then, Lemma 4.1 follows from the following. 

Lemma 4.5. Let A > 0, </> > 1 and T^ = k^. Then P^ia"" > T^) ^ as 

K OO. 

Proof. The proof relies on a restart argument. The basic idea is that 
each time the cluster at site z gets infected (between time and time T^), 
the infection dies out quickly or lives at least 2Tfi units of time. In particular, 
the number of trials before having a strong infection can be made arbitrar- 
ily large by choosing k large enough. Since a geometrical number of trials 
suffices, the first strong infection will appear by time with probability 
close to 1. To make this argument precise, we introduce, for any i>l, the 
stopping times 

Pi =inf{t > :^i(2) = 1} and pi = M{t > pi : (,tiz) = 0} 

with the convention po = 0. That is, pi is the ith time the cluster at site z 
gets infected, and pi is the ith time the cluster recovers. Let 

N = mf{i>l:pi- p,>2T^} 

so that Pat = , the first time a strong infection is born at site z. By 
decomposing the event to be estimated according to the number of trials 
N, we get 

oo 

P^ia' >T^) = J2 PdPn > T^\N = n)P^{N = n) 

n=l 

oo 

< ^ P^(pj+i — pi> Ti^/n and Pi — Pi < 2Ti^ for some i <n — 1) 

n=l 

X P^{N = n) 

oo n~l 

< X! X! ^dPi+i - Pi> ^njn and pi- pi< 2T^,)P^{N = n) 

n=l i=l 

oo 

< ^ nPg(p2 - pi > Tjn and pi - pi < 2T^)P^{N = n). 

n=l 
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To prove that the previous series can be made arbitrarily small by taking k 
large enough, we first observe that, since Xt has a drift to the right, there is 
a constant q = q{(j)) > such that 

P\{^lK) > q for any k > 1. 

This, together with Lemma 4.2, implies that -P^(to > 2Tk) = qK> 9/2 for all 
K sufficiently large. In particular, for any e > 0, there exists a large enough 
Ue such that 

oo oo 

^ nP^{N = n)=Y, nq^{l - q^^-^ < e/2. 

n=ne n=nir 

To show that the first terms tend to as k — > oo, we first observe that 

P^{P2 - Pi> T^,/n and pi - pi < 2T^,) 

< P^{p2 -pi> Tj2n) + P^{Tj2n < pi - pi < 2T^). 

Since there is at least one infected individual at site 0, we get 

PdP'i - Pi> 7)c/2n) < exp{-XT^,/2ns) for n<ne. 

In other respects. Lemma 4.4 implies that P^{T^/2n < pi — pi < 2T^) tends 
to as K — > oo. In conclusion, 

nP^{p2 - pi> T^/n and pi - pi < 2T^,) < e/2 

n=l 

for K sufficiently large, which completes the proof. □ 

Lemma 4.5 implies that if there is an infected individual at site at 
any time t G [0, T^] , then, with probability close to 1 when k is large, the 
clusters at sites z, z ~ 0, have each at least one infected individual at any 
time t € [Tf^,2Tfi], which proves Lemma 4.1. The existence of a nontrivial 
stationary distribution can then be deduced from Lemma 4.1 as in Section 3. 
□ 

5. Proof of Theorem 6. In this section we now deal with the CRP, that 
is, the recovery mechanism is now described by the transitions i — > 0, i = 
1,2,...,K, at rate 1. We prove that, contrary to the IRP, for any within 
infection rate > 0, we can find a critical value Xc{4>) such that, provided 
^ < ^ci4>), there is no epidemic whatever the cluster size. As explained in 
Section 1, the key idea of the proof is that the state of each cluster returns to 
sufficiently often so that, between two consecutive recoveries, the number 
of infected individuals in a given cluster cannot exceed some threshold kq 
with probability close to 1. 
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As previously, we rely on a rescaling argument to compare the particle 
system with the oriented percolation process introduced in Section 1. For 
any {z,n) £ Q = 'Z'^ x let Sz,n be the vertical segment and let Cz^n be 
the space-time cylinder given by (2) and (3), respectively. The good event, 
denoted by Hz^n, we now consider has to be measurable with respect to the 
graphical representation in Cz^n and assure us that, at time (n + 1)T, the 
cluster at site z is healthy. 

To construct the good event Hzfi, we first require the clusters at sites 
X, X ~ z, to recover at least once between time and time T/2 (which 
corresponds to at least one occurrence of Poisson processes with rate 1 by 
time T/2), and the cluster at site z to recover at least once between time 
T/2 and time T. The probability of this event can be estimated 

(there is at least one Xi on the segment {x} x [0,T/2] for any x ~ z 

and at least one Xi on the segment {z} x [T/2,T]) 

>l-(2(i + l)exp(-r/2), 

where Pr/ denotes the law of the GRP. In particular, given e > 0, there is a 
large enough T, fixed from now on, so that the previous event has probability 
at least 1 — e/3. The reason why we want the neighboring clusters of z 
to recover at least once by time T/2 is to control the number of infected 
individuals around z until time T; the aim is to prevent infections at site z 
coming from the outside. 

To estimate the number of infected individuals in the neighborhood of z 
until time T, we introduce the continuous-time random walk Zt G that 
makes transitions 




i + l, at rate i4>, 
1, at rate 1. 



By monotonicity of the sequence {Zt <n,yt< T}, n > 1, the Beppo-Levi 
theorem implies that 

Imi^P^iZt < n for all t<T)= P^{Zt < oo for all t<T) = l, 

where P^ denotes the law of the process Zt starting from Zq = 1. In par- 
ticular, we can find a sufficiently large kq < oo depending on (p and T, fixed 
from now on, such that 

Pzi^t > for some < i < T) < e/6d. 

Moreover, since the neighboring clusters of z recover at least once by time 
T/2, for any x ~ z and any time t G [T/2, T] , we get rit{x) < Zt in distribution 
so that 

Prj{r]t{x) > kq for some x ~ z and some T/2 <t<T)< e/3. 
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We now fix Xc{4>) so that, with high probabihty, no epidemic may occur 
provided A < Xc{4>)- When k, < kq, the result follows from Theorem 2 by 
taking Ac(</>) = Ac/kq- To deal with the nontrivial case k> kq, we complete 
the construction of our good event H^fi by requiring no infection of site z 
coming from the outside, that is, from a site x with x ^ z, until time T. Since 
each of the neighboring clusters of z has at most kq infected individuals, this 
occurs if no i-arrow with 1 <i < kq points at site z between time T/2 and 
time T, an event with probability 

exp{—d\KoT), 

which can be made greater than 1 — e/3 by choosing A > sufficiently small. 
Since there is no outside infection at site z and the cluster at z recovers 
at least once between times T/2 and T, the event assures us that z is 
healthy at time T. Finally, the probability of H^fi being greater than 1 — e 
with e > arbitrarily small. Theorem 6 follows from a comparison with the 
oriented percolation process introduced in Section 2. 
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